Spontaneous nucleation and dynamics of kink defects in zigzag arrays of trapped ions 
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The spontaneous nucleation and dynamics of topological kink defects have been studied in trapped 
arrays of 41-43 Yb ions. The number of kinks formed as a function of quench rate across the 
linear-zigzag transition is measured in the under-damped regime of the inhomogeneous Kibble-Zurek 
theory. The experimental results agree well with molecular dynamics simulations, which show how 
losses mask the intrinsic nucleation rate. Simulations indicate that doubling the ion number and 
optimization of laser cooling can help reduce the effect of losses. A range of kink dynamics is 
observed including configural change, motion and lifetime, and behavioral sensitivity to ion number. 

PACS numbers: 37.10.Rs, 52.27. Jt, 64.60.an, 64.60.Q- 
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The nucleation of topological defects during a 
symmetry-breaking phase transition has many physi- 
cal realizations. For continuous phase transitions, the 
Kibble-Zurek mechanism (KZM) provides an intuitive 
model of domain and defect formation [l], HJ and predicts 
a power-law scaling for the number of defects formed as 
a function of transition quench rate 0^3 • Several exper- 
imental tests of the KZM have been done with different 
degrees of success 0- 21 1 , all with the goal of examining 
universal behavior across diverse physical systems. 



Following a recent proposal |22j, we report on the nu- 
cleation and dynamics of topological defects in arrays of 
trapped ions, following a quench of the linear-zigzag tran- 
sition. Such studies are well suited to this experimental 
system, since it is simple enough to make direct compar- 
ison to theory, highly controllable, and efficient in that 
the defects are imaged in situ and the trapped sample 
can be recycled repeatedly across the transition. 

Laser-cooled ions held in a linear Paul trap [lH 
with strong transverse confinement organize into a one- 
dimensional (ID) linear crystal with an inhomogeneous 
axial density [Fig.QJa)]. If the transverse confinement is 
relaxed slowly, the linear ion crystal undergoes a continu- 
ous, structural phase transition to a 2D zigzag configura- 
tion, beginning at the center of the crystal [231-28] . Two 
broken-symmetry states are possible [Fig. (Ub) or (c)]. 
Rapidly quenching the transverse confinement across the 
transition can lead to crystal structures containing spon- 
taneously nucleated topological defects (kinks), which 
are formed by the interface between domains of oppo- 
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site symmetry [Figs.HKd)-(g)] (HHl, 
al. [22I, [3l[ have applied the KZM formalism to the non- 
equilibrium dynamics of the quench and calculated the 
number of kinks formed, which depends on the quench 
rate and the damping rate (77) of the system, controlled 
by Doppler laser cooling. In the under-damped regime, 
where the system's response time is dominated by the in- 
ertial dynamics of the soft transverse zigzag mode 0, [28[ , 
the rate of defect formation is independent of rj and 
given by a power law, (2<5 /2tq) & , where 2tq is the to- 
tal quench time, and 28q defines the magnitude of the 
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FIG. 1. Images of crystals of 42 174 Yb ions in various configu- 
rations: (a) ID linear crystal, (b) "zig", (c) "zag", (d) zigzag 
with a single discrete kink (indicated by arrow), (e) two dis- 
crete kinks (kink and anti-kink), and (f)-(g) single extended 
kinks. The axial trap frequency for all images is u z /2-7r = 37.6 
kHz, while the weak transverse frequency, cj v /2it, is (a) 658 
kHz, (b)-(d) 414 kHz, (e) 380 kHz and (f) 342 kHz. The 
corresponding quench end-points, Vf, are (b)-(d) 4.8 V, (e) 
5.3 V, and (f) 5.8 V. The imaging system looks at 45° to the 
transverse j/-axis where the zigzag displacement occurs. 



quench, which is linear in terms of quadratic transverse 
trap frequency. The power-law scaling is affected by the 
system inhomogeneity [22j, |32j , an effect which is largely 
unexplored experimentally. Due to the inhomogeneous 
density, the propagation speed of the transition front de- 
creases from the crystal center to the edges. This lim- 
its the size of the region where non-adiabatic dynamics, 
and so nucleation, is possible, and for the under-damped 
regime (r? 3 <C So /to) , gives a power-law scaling of b = 4/3 
(IKZM regim e) [221 . Following calculations in other con- 
texts [l5|, I34], |35|. it has been pointed out that for small 
numbers of ions or slow quenches, where the size of the 
nucleation region is smaller than the correlation length, 
the scaling doubles to b = 8/3 (DIKZM regime) 

Simultaneous experiments to ours have recently re- 
ported kink nucleation in ion arrays due to a temper- 
ature quench of a plasma cloud 37[ , and a measurement 
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of the DIKZM power law [36(. Compared to the latter 
results, our measurements of nucleation are also due to a 
quench across the structural transition but use a different 
experimental approach, including more ions and a shal- 
lower quench into the zigzag phase so that kinks retain 
their localized, discrete form. Moreover, we explore post- 
nucleation dynamics of discrete kinks including transfor- 
mation of shape, motion and lifetime in the crystal, and 
sensitivity of behavior to unit change in ion number. 

The experimental set-up is similar to that previously 
reported [38|,[39j]. We use a linear Paul trap, operating at 
a radio frequency (rf ) of fl r / /27r = 12 MHz, and composed 
of four rods with an ion-to-rod spacing of 0.7 mm and 
two end-caps separated by 2.5 mm, to form a nominally 
harmonic confining potential for ions. Photo-ionization 
is used to load a pure crystal of 42-43 174 Yb ions [39ll40j. 
which arrange as a linear string in a trap with an axial 
secular frequency of uj z /2tt — 37.6 kHz and transverse 
frequencies of u) x /2tt — 679 kHz and ui y /2ir = 658 kHz 
[Fig. Ufa)]. A single laser beam at 369.5 nm (with axial 
FWHM of 280 fiia) is used for fluorescence imaging of 
the ions and for Doppler cooling along all principal axes 
of the trap. The 369.5-nm laser is linearly polarized at 6° 
with respect to a 5.9-G magnetic field. A relatively large 
laser detuning of —33 MHz (1.5 linewidths) is used to 
suppress the effects of micromotion on fluorescence and 
cooling 4l|, arising from a non-zero component of the 
trap rf electric field along the axial direction (up to 2500 
V/m at the edges of a 43-ion crystal). With an effective 
saturation parameter of sq ~ 2.7, the estimated cooling 
rate is uniform to a factor of 2 or better across the crystal 
with variation primarily at the edges, and the calculated 
cooling time, r; -1 , at the center of the crystal is 200-400 
(is 42[ . The cooling rate places the experiment deeply in 



the under-damped regime for quench rates considered. 

To quench the system through the linear-zigzag tran- 
sition, we apply a transverse dc quadrupolar potential 
using the trap rods. The secular frequency along the 
transverse y-axis, obtained from single-ion calibrations, 
decreases with the applied dc voltage, V(> 0), according 
to uj y (t) — uj y o\/l — V/Vo where La y o/2n — 658(1) kHz 
and Vb=7.95(2) V with errors including drift. The x-axis 
frequency increases in an equivalent way. The quench 
waveform is linear (V(t) — Vf ■ £/(2tq)) with a quench 
end-point of Vf =4.2 — 5.8 V, a ramp time of 2tq = 10 — 70 
/is, and a maximum distortion of 50 mV p fc after filtering. 
The quench also incorporates small adjustments to the 
end-caps and rods to maintain a constant axial secular 
frequency and to maintain the center of the crystal at the 
rf null in all dimensions. The linear form of our quench in 
terms of Uy is as assumed in the KZM theory [22j, with 
the one modification that we use an asymmetric ramp 
across the critical point. For example, in a 42-ion string 
with w,z/27r = 37.6 kHz, where the critical transverse fre- 
quency at the center of the crystal is wjj(0)/27r = 619(3) 
kHz (or equivalently, 0.92(6) V) [43}, our quench starts at 



658 kHz (Vi =0.0 V), which ensures a sufficiently linear 
configuration [44| , and continues further past the critical 
point, for example to 414 kHz (Vf — 4.8 V) to involve 
about half of the crystal in the zigzag phase. 

Each experimental run begins with 53 ms of cooling 
in the linear configuration, followed by the quench (with 
cooling always active). Images of the initial linear and 
subsequent zigzag configurations are recorded with an ex- 
posure time of 3 ms. A 200- /is ramp recycles the crystal 
to the linear structure. Data is acquired in batches of 
100-700 runs over 3 minutes. Our data sets contain 80% 
of runs approximately equally divided between 42 and 43 
ions with 41 or fewer ions in the remainder. Occasionally, 
the fluorescence of an ion is interrupted by internal-state 
pumping or molecule formation through background-gas 
collisions [3^] . Runs with dark ions within the zigzag re- 
gion of the crystal (< 14% of runs, typical) are excluded. 
Image analysis is used to identify locations of any kinks 
in the zigzag structure. For configurations without kinks, 
temporal statistics of the two symmetry-broken phases 
are nominally random with a bias of <1% and a normal- 
ized binary autocorrelation on adjacent runs of <1%. 

In the first experiment, we explore the formation of 
two types of kink (2t|: discrete kinks, which are localized 
to about 3 ion-sites [Figs. QJd),(e)], and extended ones, 
which involve more sites and occur in crystals brought 
far into the zigzag region [Figs. [TJf),(g)]. Figures EJa)- 
(c) show the observed axial distribution of both types and 
their total numbers as a function of the quench end-point, 
Vf. The data, for 41-43 ions, is taken at a constant ramp 
rate of 0.19 V//iS so that the nucleation dynamics near 
the critical point is expected to be same. Any differences 
are due to post-nucleation dynamics during the remain- 
der of the ramp through to detection. To understand the 
outcome of the experiment, we have performed molecular 
dynamics simulations with no free parameters: The ion 
trap is treated as an ideal harmonic potential limited to 
two dimensions. The small axial anharmonicity in the 
experiment is modeled in an effective way by shifting the 
axial frequency from 37.6 to 38.2 kHz [43|. The effect 
of the Doppler cooling laser is modeled explicitly 45 1 
and includes the steady-state fluorescence behavior of 
174 Yb + and all laser calibrations excluding beam pro- 
file |38[. Simulations for 42 ions reveal similar distribu- 
tions to Figs.[^a) and (b). There is also a good qualita- 
tive match in the numbers of discrete and extended kinks 
between simulation and experiment [see Fig. [21(c)]. From 
the simulations, we find that discrete kinks are nucleated 
in the central 5-10 sites of the crystal. Then, there is an 
outward motion and loss of the discrete kinks, propelled 
by the gradient in zigzag amplitude and axial density, 
which affect the underlying potential for the kinks. As 
the quench end-point increases towards 5.0 V, the num- 
ber of discrete kinks and the width of their distribution 
increase since more of the dispersing kinks are trapped 
in a deeper and wider zigzag region. Above 5.0 V, the 
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FIG. 2. (color online) Observed axial position distribution of 
(a) discrete and (b) extended kinks versus quench end-point. 
All data taken at a quench rate of 0.19 V//zs. (c) Average 
number per run of discrete kinks (red circles) and extended 
kinks (blue squares) for distributions in (a) and (b). Data 
points, for 41-43 ions, include 1000-2200 runs. Errors calcu- 
lated assuming binomial statistics. Lines are a simulation for 
42 ions with no free parameters, and shown 2 ms after the 
quench, 2/3 of the camera exposure time. Statistical uncer- 
tainties in simulation comparable to experiment ones shown, 
(d) Simulation frames showing evolution after nucleation of a 
discrete kink (red) into an extended kink (blue) at Vf = 5.8 V 
(oi !/ /27r=342 kHz). Middle 144 fim of a 42-ion crystal shown. 



number of discrete kinks collapses when, as seen in the 
simulation, those at the center of the crystal distort to 
extended form. Extended kinks can also migrate to the 
center [Fig. HJd)], suggesting a change in the underlying 
potential for the kinks and leading to the well-centered 
distribution observed in the experiment [Fig. [2jb)]. Due 
to reduced losses, the number of extended kinks in the 
simulation saturates at a value 1.4 times higher than the 
peak in discrete number. The experiment data ends at 
5.8 V, well before the ~7.1-V onset of another structural 
transition to a different 2D configuration. As this tran- 
sition is approached, we see an enhancement in (other- 
wise negligible) background kink nucleation, due to en- 
ergy transfer events such as background-gas collisions. 

Next, we proceed to study kink nucleation as a func- 
tion of quench rate. With an eye towards subsequent 
dynamics studies, we work with discrete kinks and op- 
erate at Vf = 4.8 V, where their maximum number is 
detected [Fig. [2jc)]. Figure [3ja) shows a log-log plot of 
the average number of kinks detected per run (black tri- 
angles) as a function of scaled rate, 1/(uj z ■ 2tq). Data 
for 41-43 ions is combined. As the quench rate increases 
from 1/(cl>2 • 70/is), the number of kinks grows, peaks 
at l/(id z • 15 fJ.s) and begins to drop. Figure [3ta) also 
includes results of molecular dynamics simulations, per- 
formed as described above. The first simulation (solid 
black line) shows the number of kinks 2 ms after the 
quench (2/3 of the exposure time in the experiment), by 
which time the kink number has stabilized. The simu- 
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FIG. 3. (color online) (a)Average number of kinks observed 
per run vs. scaled quench rate for standard trap with ui z /2ti = 
37.6 kHz (triangles), and for one with secular frequencies re- 
duced by 2 (red squares) and quench changed accordingly. 
Data, for 41-43 ions, include 1000-5000 runs. Error bars 
are calculated binomial, approximately consistent with fluc- 
tuations. Solid line is a molecular dynamics simulation with 
parameters including uj z /2tv — 38.2 kHz (see text), and shown 
2 ms after the quench. Statistical uncertainties of simulation 
less than or equal to the experiment. Dashed black line is the 
same simulation with kinks counted shortly after nucleation 
(>7 fis). To suppress thermal fluctuations, the kink search is 
limited to the zigzag region with minimum amplitude of 0.54 
/»m. Dotted and dash-dotted lines (green) show kink number 
2 ms after quench and shortly after nucleation, respectively, 
for a simulation of 84 ions with optimized Doppler cooling 
(detuning -11.3 MHz, s = 2.0) and with uj z /2tt = 20.6 kHz, 
chosen to match the transverse critical frequency for 42 ions 
at crystal center. Simulation plot scaled using uj z /2-k = 38.2 
kHz for direct comparison to 42 ions, (b)-(d) Observed po- 
sition distributions of kinks for 42 ions (red, upper) and 43 
(blue, lower) in standard trap at quench time 2r 9 of (b) 20, 
(c) 15 and (d) 12 ^is with simulations (lines) for comparison. 



lation matches the observed number of kinks extremely 
well, and reproduces features in the observed kink dis- 
tributions, related to rapid, post-nucleation dynamics 
[Fig.[3jb)-(d)]. To understand the effect of kink losses be- 
fore detection in the experiment, we use the simulation to 
assess the "early-counted" number of kinks shortly after 
nucleation (dashed black line). At fast quenches, there 
is enhanced loss due to kinetic energy from nonadiabatic 
excitation of vibrational modes, and before the crystal is 
re-cooled. Kink pair annihilation also plays a minor role. 
At slow quenches, we again see a loss of kinks, which 
worsens with decreasing quench rate. In this case, the 
zigzag region expands too slowly to trap the kinks, and 
they escape to the zigzag edges during the quench. 

Although kink losses in the experiment mask the in- 
trinsic nucleation rate - see, however, Ref. (36| - we com- 
ment briefly on the power law scaling for the simulation 
and experiment. Figure |3ta) covers about three octaves 
in quench rate separated by log 2 (l/2o; 2 TQ) = — 3.1 and 
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2.1. The KZ theory [22J, |36J suggests that our data 



covers two regimes of nucleation: For the slow-quench 
range, we expect a DIKZM power law of 8/3 = 2.67, which 
agrees with the fit value of 2.5(2) for the early-counted 
simulation (with significant systematic error associated 
with rapid losses). For the fast-quench range, where the 
nucleation region is large enough that nucleation of two 
or more kinks becomes significant in the simulation, one 
expects to cross over towards the IKZM power law of 
4/3 = 1.33. The early-counted simulation gives a fit value 
of 1.7(1). We note that all cases of two kinks observed 
in the experiment [Fig. [lid)], albeit rare, occur in or ad- 
jacent to the fast-quench region. Finally, in the middle 
octave of quench rate, where the kink losses are least, a 
fit to the experimental data gives a power law of 3.3(2), 
compared to the early-counted simulation value of 2.0(1). 

For comparison, we also measure kink nucleation in a 
trap with all secular frequencies reduced by a factor of 2 
(red squares in Fig.[3f a)). Overall, the data match results 
at higher trap frequency, which shows the weak sensitiv- 
ity of nucleation and losses to effective cooling rate, rj/u} z , 
over such a variation, except perhaps for losses at the 
fastest quenches [46j. Of several possibilities to reduce 
the losses in our experiment, we note that an increase in 
ion number can help, particularly for slow quenches by 
reducing the crystal inhomogenieties that transport dis- 
crete kinks outward. As shown in Fig.^a), a doubling of 
ion number (and using twice better cooling) already leads 
to significant reduction in losses over the same quench 
range, albeit in a scaling regime closer to IKZM. 

We conclude with an investigation of the dynamics 
of single discrete kinks, and measure their axial motion 
in the crystal and their lifetime. The experiments are 
performed at V) = 4.8 V and 2tq = 20 /is, for which 
the initial distribution of the kinks is relatively uniform 
[Fig. EJb)]. Figures Ufa) and (b) show the trajectories of 
kinks in crystals of 42 and 43 ions for three initial axial 
positions. The kinks move between sites corresponding to 
local potential wells separated by the axial lattice spac- 
ing (5.6 /zm) [3(|. The 40-ms frame rate of our camera 
does not permit us to resolve fast dynamics, including 
the motion between sites. Nevertheless, it can be seen 
that, for 42 ions, most of the kinks that begin at the cen- 
ter of the crystal either stay there or move by one site. 
In comparison, the central kinks in a string of 43 ions 
migrate to the edges faster. For both cases, kinks closer 
to the edges tend to move outwards, as expected. 

In Figs. |4|c) and (d), we show the exponential life- 
time of kinks as a function of their initial position as 
well as the distributions of kink positions at t = s and 
t = 0.4 s. In Fig. Etc), we can see that the lifetime of 
kinks for 42 ions is highest near the middle of the crys- 
tal (with the ~2-s value likely limited by background-gas 
collisions) , and decreases roughly symmetrically towards 
the edges. Losses due to the position-dependent lifetime 
largely account for the observed time evolution of the po- 
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FIG. 4. (color online) (a)-(b) Axial trajectories of discrete 
kinks in zigzag crystals (V/ = 4.8 V) with (a) 42 and (b) 43 
ions for three different initial positions: center and four sites 
away on both sides. The number of trajectories for each case 
is 11 with fewer for z < 0. (c)-(d) Exponential lifetime of 
kinks (filled circles) vs initial position in crystals of (c) 42 
and (d) 43 ions. Lifetimes and their uncertainties obtained 
from a maximum likelihood estimate including binning due to 
camera frame rate. Also shown are axial position distributions 
at t — ms (filled gray bars) and t = 400 ms (open red bars). 



sition distribution. In Fig. HJd) , for 43 ions, the 0.3-0.5 s 
lifetime of central kinks contrasts with their complete ab- 
sence in the position distribution after 0.4 s, and points 
to their migration away from the deeper zigzag region at 
the center of the crystal [47[ . Simulations show that this 
is associated with new dynamics in which the kink hops 
two sites (as observed in Fig. 0|b) center) via a transient 
extended-like form. The difference between 42 and 43 
ions is understandable, considering that the experiment 
is performed at Vf — 4.8 V, near to the sharp onset of the 
bimodal distribution in Fig. [2fa). Since the zigzag phase 
begins 0.3 V lower for 43 ions, the onset of the bimodal 
distribution will also occur at a lower voltage. 

To conclude, we have measured the number of kinks 
formed in ion crystals as a function of quench rate in the 
under-damped regime of the KZM, and have shown how 
losses mask the intrinsic nucleation rate. Although simul- 
taneous results based on a deeper quench and extended 
kinks show lower loss [36[, our method of nucleation as- 
sessment should be advantageous for larger ion num- 
bers with multiple kinks, since the localized character of 
the discrete kinks will help reduce interaction-dependent 
losses. We have also studied the post-nucleation dy- 
namics of kinks, and have observed long lifetimes, which 
may allow future explorations of internal kink dynam- 
ics [29|, [30[. A natural future direction is the study of 
linear-zigzag dynamics in the quantum regime [48j . 
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